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Abstract 

Within the framework of generalized factorization of higher-twist contributions, including modification to splitting 
functions of both quark and gluon, we get and numerically resolve the medium-modified DGLAP (mDGLAP) evolu- 
tion equations. With Woods-Saxon nuclear geometry and Hirano 3D ideal hydrodynamic simulations of hot medium, 
we study the medium modified fragmentation functions (mFF) in DIS and Au+Au collisions in RHIC. Our calculation 
imply that the parton density in hot medium produced in RHIC is about 30 times larger than cold nucleon. 
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1. Introduction 

Since jet quenching phenomenons are observed in RHIC, many phenomenological studies of it indicate a scenario 
of strong interaction between energetic partons and the hot medium with an extremely high initial parton density 
The same phenomena are also predicted in deeply inelastic scattering (DIS) off large nuclei when the struck quark 
propagates through the target nuclei J3]. 

In the presence of nuclear or hot QCD medium, the initially produced energetic partons will have to go through 
multiple scattering and induced gluon bremsstrahlung before hadronization. The induced gluon bremsstrahlung ef- 
fectively reduces the leading parton's energy and softens the final hadron spectra or parton fragmentation functions. 
To take into account multiple induced gluon emissions, one can follow the resummation of gluon bremsstrahlung in 
vacuum and assume that multiple medium induced bremsstrahlung can be resummed in the same way to obtain the 
mDGLAP evolution equations for the modified fragmentation functions llllz 131. 
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where the modified splitting functions are given by the sum of the vacuum ones and the medium modification 
l\) = Ja^bdz) + Ay a _» fec (z, fij), 
The first mDGALP evolution equation for quark in Eq. ([T) is derived in Refs. (H 
modification to the splitting functions for quark 



Ml in DIS off nuclei. The 
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are obtained from the induced gluon bremsstrahlung spectra and therefore are related to the twist-four nuclear quark- 
gluon correlation distribution which essentially describe the amplitude of the second hard parton scattering and the 



induced gluon radiation. The matrix element A^- q 'T qg in the second term in Eq.© comes from virtual corrections to the 
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multiple parton scattering cross section. This term can be constructed from the momentum sum rule (or momentum 
conservation) for mFF l!l|2|,|5|], J dzzAD h q (z, Q 2 ) = 0. 

In order to get the complete coupled mDGLAP evolution equations, we also consider multiple scattering and 
induced gluon bremsstrahlung for a gluon jet to get the mDGLAP evolution equation for gluon in Eq. (0. From 
gluon-gluon scattering matrix elements, one can obtain the medium modification to the splitting functions for gluon 



Ajg^ g g and Ay g ^ qq (see details in Ref. 1 14]). 

In terms of the generalized jet transport parameter q in nuclear medium, we can express the quark-gluon correlation 
function T q approximately as the integration over the parton's trajectory through the medium 
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dy-q F (y)4 sm 2 (x L p + y- /2) . (5) 



Given the initial conditions of fragmentation functions at initial energy scale D a {Qfy, We can numerically solve 
the coupled mDGLAP evolution equations using modified HOPPET I Kill Fortran 95 package in LO. Such initial 
conditions in principle should be different in medium and vacuum. To take into account medium modification to 
the fragmentation functions at the initial scale Qq, we will assume in this study fl3ll ZX.fO 2 ,) = D a (Qfy + ADaiQfy, 
where D a {Qfy is the vacuum fragmentation function and AD a (Qfy is generated purely from medium via the mDGLAP 
starting at p 2 = 0. 



2. Modified Fragmentation Function in DIS 



To calculate the mFF in semi-inclusive DIS off a nucleus, we employ the Woods-Saxon nuclear geometry. Consid- 
ering the initial quark jet produced at yo that travels along a direction with impact parameter b, we assume that the jet 
transport parameter along the quark jet trajectory is proportional to the nuclear density q(y, b) = qopAiy, b)/pA(0, 0). 

If we neglect the nuclear and impact parameter dependence of the nuclear quark distribution function, the photon- 
nucleon cross section that produces a quark at (yo, b) is proportional to the nuclear density distribution pAiyo, b). Then 
the averaged mFF should be 
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In order to calculate the D(z,yo, b) for a quark produced at location (yo, b), the path integral in the modified splitting 
functions should be replaced by the following 
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Shown in Fig. (Q} are the calculated nuclear modification factors of n ± , and p(p) with different values of jet 
transport parameter qo = 0.015 + 0.005 GeV 2 as compared to the HERMES experimental data fTlll : 



fljffcv) 



N h (z,v) 
N e (v) 



u / 



N h (z,v) 
N e (v) 



(Ze 2 f q f (x)D h f (z) 



Ee 2 f q f (x) 



(Ze 2 f q f (x)D h f (z) 



L 

Y.e 2 . 



jq f (x) 



(8) 



Illustrated in Fig. (Q]), the medium modification of the mFF gradually disappears as the initial jet energy E in- 
creases. The agreement between our theory calculations and experimental data is generally good except at lower 
energy where hadronic absorption might become important. We can see significant different modification at HER- 
MES for p and p, this may be because of the process of quark-antiquark annihilation in twist-4 double scattering and 
the asymmetry of q and q in nuclei C1QJ], which we have not put into mDGLAP yet. 

As we have discussed, the initial condition for mFF in the medium at = 1 GeV 2 is different from which in 
vacuum. Therefore, most of the medium modification to the mFF come from mDGLAP evolution at low Q 2 while 
contribution from high Q 2 region is power-suppressed. This will lead to a very weak Q 2 dependence as shown in the 
left panel of Fig. [2] The calculated suppression factors are almost independent of Q 2 , consistent with the experimental 
data 111 ill . If one has chosen the initial condition at Qq as the same as the vacuum one, one would obtain a too strong 
Q 2 dependence of modification factor. 
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3. Modified Fragmentation Function in Au+Au Collisions 

One can extend the calculation for medium modified parton fragmentation functions in DIS to hot medium 
like QGP or hot hadronic matter created in high-energy heavy-ion collisions. To take into account both the the 



longitudinal and transverse expansion of the hot matter we use a 3D ideal hydrodynamic simulations 11181, 119H which 
give us information of the temprature, energy density, the fraction of the hadron phase and so on, on each step of the 
hot matter evolution. 

The jet transport parameter in hot medium at given time t and local position r is assumed to include the contribu- 
tion from both the QGP phase and hadronic phase 1 1 7| 

M . . P QCP (T,r) P/ ,(r,r) 

4(r,r) = 4o-^7 7zO--f) + <lN /. (9) 

P yw TTo,0) p N 

where f is the fraction of the hadronic phase, qo is the jet transport parameter at the center of the bulk medium in the 
QGP phase at its formation time tq. p/,(r, r) is the number density of the hadron resonance gas, p# = no ~ 0.17 fm~ 3 
is the nucleon density in the center of a large nucleus and = 0.015 GeV 2 /fm is the jet transport parameter in cold 
nuclear matter we got in last section. 

Then we can calculate the mFF D h p {zh, Q 2 , E, r, <p, b) in hot medium, and average it over the initial parton produc- 
tion position and the out going directions. Assume the parton production cross section is proportional to the overlap 
function, we can get 
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If we neglect nuclear effect such as the shadowing effect on the initial parton distribution function, we can get the 
nuclear modification factor Raa for a fixed impact parameter b 

_ dcr' AB /dyd 2 p T d 2 b _ fPQci , x 2 ) g do- ® (gfa, Q 2 , E, b)> 
RaaO - TAA[b)d(r h p/dyd 2 pT ~ f P ( XuX2 )®dcT®D(z h ,Q 2 ) ' (U) 

Shown in the right panel of Fig. [2] is the comparison of our result about the medium modification factor for top 5% 
centrality Au+Au collisions to the PHENIX data [20]. The jet transport parameter at the center of the bulk medium 
in the QGP phase at its formation time qo = 0.5 + 0.1 GeV 2 /fm can fit the experiment data in error bar. 
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Figure 2: left: Comparison of the modified multiplicity ratios as a function of Q 2 at fixed value of z and jet energy E with the HERMES fllll data 
for Ne, Kr and Xe targets. For clear presentation the modification factors for different targets have been shifted vertically by some value(A> by 
-0.05 and Xe by -0.2). right: comparison of the nucleon modification factor Raa for n° in - 5% biased events. 



4. Conclusions 

We have got the coupled mDGLAP evolution equations within the framework of generalized factorization of 
higher-twist contributions. By solving the mDGLAP equations numerically, we study the nuclear modified fragmen- 
tation functions both in cold nucleon and hot medium produced in RHIC. From our results, one is suggested that the 
gluon density in QGP is about 30 times larger than which in cold nucleon. 
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